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ABSTRACT 

Over the duration of the Kepler mission, KIC 8462852 was observed to undergo irregularly 
shaped, aperiodic dips in flux of up to ^ 20%. The dipping activity can last for between 5 
and 80 days. We characterize the object with high-resolution spectroscopy, spectral energy 
distribution fitting, radial velocity measurements, high-resolution imaging, and Fourier anal¬ 
yses of the Kepler light curve. We determine that KIC 8462852 is a typical main-sequence 
F3 V star that exhibits no significant IR excess, and has no very close interacting compan¬ 
ions. In this paper, we describe various scenarios to explain the dipping events observed in 
the Kepler light curve. We confirm that the dipping signals in the data are not caused by 
any instrumental or data processing artifact, and thus are astrophysical in origin. We con¬ 
struct scenario-independent constraints on the size and location of a body in the system that 
is needed to reproduce the observations. We deliberate over several assorted stellar and cir- 
cumstellar astrophysical scenarios, most of which have problems explaining the data in hand. 
By considering the observational constraints on dust clumps in orbit around a normal main- 
sequence star, we conclude that the scenario most consistent with the data in hand is the 
passage of a family of exocomet or planetesimal fragments, all of which are associated with a 
single previous break-up event, possibly caused by tidal disruption or thermal processing. The 
minimum total mass associated with these fragments likely exceeds 10“® Mq, correspond¬ 
ing to an original rocky body of > 100 km in diameter. We discuss the necessity of future 
observations to help interpret the system. 

Key words: stars: individual (KIC 8462852), stars: peculiar, stars: activity, comets: general, 
planets and satellites: dynamical evolution and stability 


1 INTRODUCTION 

For over four years, NASA’s Kepler mission measured the bright¬ 
ness of objects within a ~ 100 square-degree patch of sky in the 
direction of the constellations Cygnus and Lyrae. The program’s 
targets were primarily selected to address the Kepler mission goals 
of discovering Earth-like planets orbiting other stars. Kepler tar¬ 
geted over > 150,000 stars, primarily with a 30-minute observing 
cadence, leading to over 2.5-billion data points per year (> 10 bil¬ 
lion data points over the nominal mission lifetime). 

The Kepler mission’s data processing and identification of 
transiting planet candidates was done in an automated manner 
through sophisticated computer algorithms (e.g., Jenkins et al. 
2010). Complementary to this analysis, the Zooniverse citizen sci¬ 
ence network provided the means to crowd source the review of 
light curves with the Planet Hunters project^ (e.g., Fischer et al. 
2012). In this framework. Planet Hunter volunteers view 30 day 
segments of light curves in the ‘Classify’ web interface. A vol¬ 
unteer’s main task is to identify signals of transiting planets by 
harnessing the human eye’s unique ability for pattern recognition. 
This process has shown to have a detection efficiency to identify 
planetary transits > 85% using the first Quarter of Kepler data 
(Schwamb et al. 2012). The Planet Hunters project has now discov¬ 
ered almost a hundred exoplanet candidates, including several con- 


* Based on observations obtained with the Nordic Optical Telescope, oper¬ 
ated on the island of La Palma jointly by Denmark, Finland, Iceland, Nor¬ 
way, and Sweden, in the Spanish Observatorio del Roque de los Muchachos 
of the Instituto de Astrofisica de Canarias. 

f The data presented herein were obtained at the W.M. Keck Observatory, 
which is operated as a scientific partnership among the California Institute 
of Technology, the University of California, and the National Aeronautics 
and Space Administration. The Observatory was made possible by the gen¬ 
erous financial support of the W.M. Keck Foundation. 

1 www.planethunters.org 


firmed systems (Fischer et al. 2012; Lintott et al. 2013; Schwamb 
et al. 2013; Wang et al. 2013; Schmitt et al. 2014). 

Because Planet Hunter volunteers look at every light curve 
by eye, serendipitous discoveries are inevitable, especially in rich 
data sets such as that which Kepler has provided. As such, a key 
aspect of the Planet Hunters project is the ‘Talk’ interface. ‘Talk’ 
is a community discussion board/site where volunteers can discuss 
light curves and present further analysis on objects viewed in the 
main ‘Classify’ interface. In a handful of cases, such as the discov¬ 
ery of the unusual cataclysmic variable, KIC 9406652 (Gies et al. 
2013), the default aperture mask used to generate the Kepler light 
curve was not perfectly centered on the object of interest. Because 
of this, interesting events in the Kepler light curve would appear to 
come and go as a result of the shifting orientation of the aperture 
mask when the spacecraft underwent a quarterly rotation. Events 
such as these are tagged and discussed on ‘Talk’, making it pos¬ 
sible to return to the raw data target pixel files (TPF) to extract 
improved light curves with modified aperture masks, for example. 

This paper presents the discovery of a highly unusual dip¬ 
ping source, KIC 8462852, from the Planet Hunters project. In just 
the first quarter of Kepler data. Planet Hunter volunteers identi¬ 
fied KIC 8462852’s light curve as a “bizarre”, “interesting”, “giant 
transit” (Q1 event depth was 0.5% with a duration of 4 days). As 
new Kepler data were released in subsequent quarters, discussions 
continued on ‘Talk’ about KIC 8462852’s light curve peculiarities, 
particularly ramping up pace in the final observations of the Kepler 
mission. 

In this work we examine the full 4 years of Kepler observa¬ 
tions of KIC 8462852 as well as supplemental data provided by 
additional ground- and space-based observations. In Section 2, we 
characterize KIC 8462852 using Kepler photometry, spectroscopic 
analysis, AO imaging, and spectral energy distribution analysis. We 
discover a wide M-dwarf companion to the system and argue that 
with the data sets we have in-hand, we can exclude the presence 
of an additional massive gravitationally bound companion nearby. 
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In Section 3, we visit possible explanations for the peculiar ob¬ 
servations of KIC 8462852, including instrumental artifacts, intrin¬ 
sic/extrinsic variability, and a variety of scenarios invoking light¬ 
blocking events. We formulate a variety of scenario independent 
constraints in Section 4, and elaborate on specific occultation sce¬ 
narios in Section 5. In Section 6, we conclude by discussing future 
observations needed to conshain the nature of the object. 


2 DATA 

KIC 8462852, also known as TYC 3162-665-1 and 2MASS 
J20061546-1-4427248, is a 17 ~ 12 mag star in the Kepler field of 
view. Its light curve was identified serendipitously by the Planet 
Hunters project, and was deemed an interesting object that was 
worthy of further investigation. In the following sections, we char¬ 
acterize the system with data from Kepler as well as additional data 
from various targeted and archived programs. 



Frequency (cycles/day) 


2.1 Kepler photometry 

The Kepler mission was launched on 2009 March 7, and it started 
science observations on 2009 May 13. The nominal mission was 
finished almost 4 years later, on 2013 May 12, after the failure of 
the second reaction wheel. KIC 8462852 was observed throughout 
the main Kepler mission (divided into Quarters 0 - 17) under long- 
cadence (30-minute) observations yielding an ultra-precise, nearly 
uninterrupted, light curve during this time. Kepler data files pro¬ 
vide both the ‘uncorrected’ Simple Aperture Photometry (SAP) and 
the ‘corrected’ Pre-search Data Conditioning (PDCSAP) fluxes (for 
details, see Christiansen et al. 2012). In this work, our analysis uses 
the normalized, PDCSAP data. Note that we have thoroughly val¬ 
idated the data to ensure that any flux variations represent phys¬ 
ical events in or near the star (and they do); these processes are 
described in detail within Section 3.1, and we do not repeat them 
here. 

In Figure 1, we present a montage of plots capturing much 
of the interesting flux variations observed in the Kepler timeseries 
data. The top two panels, ‘(a)’ and ‘(b)’, show the flux time series 
for the entire Kepler mission, but with different vertical flux scales. 
These show that the flux is relatively constant for most of that time, 
but is punctuated by a number of substantial dips in flux, includ¬ 
ing a 15% drop near day 800, and a whole sequence of dips (with 
one reaching a depth of 22%) after day 1500. Panel ‘(b)’ marks 
the occurrence of 10 discrete dips (see Table 1). For convenience, 
we hereafter refer to the two main dip structures between day 788 
and 795 and between day 1510 and 1570, as events ‘D800’ and 
‘D1500’, respectively. Panel ‘(c)’ is a zoom in on the dip D800. 
The remaining three panels are progressively zoomed in around the 
exotic complex of dips at D1500. 

The D800 dip feature is clean, sharp, and asymmetric in shape. 
It possesses a gradual dimming lasting almost a week, and transi¬ 
tions back to its nominal brightness in just a couple of days. The 
D1500 complex consists of many dips, with variable shape and du¬ 
ration, often occurring concurrently as if several independent oc¬ 
cultation events were superimposed upon each other. The D1500 
dips persist for ~ 100 days until the Kepler mission’s end, and 
only for a small part of this time does it appear ‘quiescent’. There 
are also other smaller ~ 0.5% dips, including three earlier in the 
mission around day 140, day 260, and day 359, and another after 
the D800 event, around day 1205 (dips #1, 2, 3 and 6, respectively; 
Figure 1 ‘(b)’. Table 1). Several more 0.5 — 1% dips appear in 


Figure 2. Fourier transform for KIC 8462852. The peaks are labeled with 
the harmonic numbers starting with 1 for the base frequency. Refer to Sec¬ 
tion 2.1 for details. 

and around the two deep D1500 features, including a ~ 3% dip 
around day 1540. Two small dips occurring at day 1205 and day 
1540 have shapes with a similar distinctive, ‘triple-dip’, symmet¬ 
ric profile, however, they differ in duration by a factor of 3 and in 
degree of dimming by a factor of 5. All of the fluctuations in in¬ 
tensity visible on these plots are real, i.e., not due to statistical or 
instrumental variations (Section 3.1). 

There are also modulations in the raw flux data at the ~ 
500 ppm level which are visible by eye. To further explore whether 
any of these modulations are periodic, or have a periodic compo¬ 
nent, we generated a Fourier transform (FT) of the data with the 
dips excised from the data train. Figure 2 shows the FT of the Ke¬ 
pler photometry and one can see a clear periodicity of 0.88 day 
(1.14 cycles/day) and its next two higher harmonics. 

This 0.88-day signal is a broad feature that resembles typical 
FTs of Kepler targets for early type stars (Balona 2013, see their 
figure 6). If this is a rotation period, then the projected rotational 
velocity (from Section 2.2) of 84 ±4 km s“^ represents a minimum 
stellar radius of ~ 1.46 R©, consistent with the radius of an F-type 
star (also see Section 2.2). Also seen in Figure 2 just to the left 
of the base frequency is a broad collection of smaller peaks. This 
suggests that something more complicated than a single rotating 
surface inhomogeneity is producing the observed signal. 

We investigate the stability of the frequencies observed in the 
FT by performing a Short-Term Fourier Transform (STFT), again 
clipping the data in the dipping regions. In the STFT method, the 
data are broken up into “short” segments of 43 d. This segment du¬ 
ration has been selected to optimize both time and frequency reso¬ 
lution. The FT is computed and displayed vertically on the plot, and 
this is repeated as a function of time, with overlap in time segments 
to gain back some temporal resolution. 

The STFT is presented in Figure 3. This shows that the 
0.88 day signal is present in most of the Kepler time series, with 
the strongest presence occurring around day 1200. Interestingly 
however, around day 400 and day 1400, we see major contribu¬ 
tions at different frequencies, corresponding to ~ 0.96 days and 
~ 0.90 days, respectively. We conclude that these are the source of 
the broad collection of peaks to the left of the base frequency noted 
above. These low-frequency side-bands could possibly be due to 
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Time (BJD-2454833) 





Figure 1. Montage of flux time series for KIC 8462852 showing different portions of the 4-year Kepler observations with different vertical scalings. The top 
two panels show the entire Kepler observation time interval. The starting time of each Kepler quarter Q is marked and labeled in red in the top panel ‘(a)’ • Dip 
numbers corresponding to the 10 discrete dips listed in Table 1 are marked and labeled in blue in panel ‘(b)’. Panel ‘(c)’ is a blowup of the dip # 5 near day 
793 (D800). The remaining three panels, ‘(d)’, ‘(e)’, and ‘(f)’, explore the dips (labeled in blue) which occur during the 90-day interval from day 1490 to day 
1580 (D1500). Refer to Section 2.1 for details. 
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Kepler time (days) 

Figure 3. The STFT for the Kepler flux time series. The main base period of 
~ 0.88 days is present throughout the span of observations. We identify (at 
least) two additional frequencies appearing around day 400 and 1400, cor¬ 
responding to periods of 0.96 to 0.90 days, which we attribute to differential 
rotation. Refer to Section 2.1 for details. 


regions contrasted in flux (e.g., starspots, chemically peculiar re¬ 
gions) appearing at higher latitudes coupled with differential ro¬ 
tation. This is consistent with the differential rotation (or inferred 
fractional frequency difference of ~ 10%) for F-type stars (Rein¬ 
hold et al. 2013). We would like to note however, that we cannot 
completely discount the possibility that these periods are due to 
pulsations. The position of KIC 8462852 is within the Gamma Do- 
radus (7 Dor) region of the instability strip, where pulsations are 
observed at < 5 cycles d“^ (e.g., Uytterhoeven et al. 2011). To 
investigate this, we then compared the STFTs of known 7 Dor pul- 
sators to the STFT of KIC 8462852. We found that the dominant 
frequencies in STFTs for known 7 Dor stars do not evolve with 
time, contrary to the STFT for KIC 8462852. This supports the in¬ 
terpretation that the ~ 0.88 d signal is due to the star’s rotational 
period. 

We also report on the presence of variability on the timescale 
of 10 - 20 days (Figure 2), which, when present, is visible by eye 
in the light curve^. We illustrate this in Figure 4, showing zoomed 
in regions of the Kepler light curve. The star’s 0.88 d period is also 
evident in each panel as the higher-frequency flux variations. The 
panel second from the bottom ‘(c)’ shows no low-frequency (10 
- 20 day) variations, but the rest do. While the largest of the dip¬ 
ping stmctures within the D1500 events could also be described as 
having a periodic stmcture close to 20 days, the magnitude of the 
variability and the temporal behavior are much different than these 
low-amplitude variations described here. Thus, we cannot suggest 
any connection between the D1500 features and the 10-20 day vari¬ 
ability. Finally, we note that the 10-20 day variability may actually 
arise on a faint neighboring star (see Sect. 2.3). 

There is another possible periodicity that is worth discussing 
briefly. In Table 1, we summarize the times and depths of 10 dis¬ 
crete dips present in the Kepler light curve, also labeled in panel 
‘(b)’ - ‘(e)’ of Figure 1. If we examine the two most prominent 
dips (DI568 and DI520; also see panel ‘(d)’ in Figure 1), we see 
that they have a separation of ~ 48.8 days. We can also see that 
the D800 dip (dip #5 in Table 1) is separated from the D1520 dip 
by 15 of these intervals, if the interval is more precisely defined to 
be 48.4 days. Furthermore, the very shallow dips early in the Ke¬ 
pler time series at D260 and D360 are very close to 26 and 24 of 
these 48.4-day cycles from the D1520 dip. The other five identified 
discrete dips (four of which are very shallow), also listed in Ta¬ 
ble 1 , are about a half cycle out of phase with this period to within 


^ Also present in the raw SAP data. 






Figure 4. Stacked plots showing a zoomed-in portion of the Kepler light 
curve. The star’s rotation period of 0.88 d is seen in each panel as the high- 
frequency modulation in flux. With the exception of panel ‘c)’, a longer 
term (10 -20 day) brightness variation is observed, also present in the FT 
shown in Figure 2. Refer to Section 2.1 for details. 


Table 1. Principal Dip Times of KIC 8462852 vs. 48.4-day Period 


dip # 

name 

depth 

BID 

(- 2 454 833) 

cycles 
(from dip 5) 

|residual| 
(from integer) 

1 

(D140) 

0 . 5 % 

140.49 

-13 

0.52 

2 

(D260) 

0 . 5 % 

261.00 

-11 

0.01 

3 

(D360) 

0 . 2 % 

359.11 

-9 

0.04 

4 

(D425) 

0 . 2 % 

426.62 

-7 

0.44 

5 

(D800) 

16% 

792.74 

0 

0.00 

6 

(D1200) 

0.4% 

1205.96 

8 

0.54 

7 

(D1500) 

0.3% 

1495.97 

14 

0.53 

8 

(D1520) 

21% 

1519.60 

15 

0.02 

9 

(D1540) 

3% 

1540.40 

15 

0.45 

10 

(D1570) 

8% 

1568.49 

16 

0.03 
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~ ±5% of a cycle. In this exercise, we have neglected the fact that 
the three most prominent dips in the D1500 region are quite highly 
structured, and they also have additional minima whose times could 
have been tagged and included in the analysis. At this time we do 
not ascribe any particular significance to this period, but it is some¬ 
thing to bear in mind as more data on this object become available. 

2.2 Spectroscopy 

We obtained four high resolution (R — 47000) spectra of 
KIC 8462852 with the FIES spectrograph (Frandsen & Lindberg 
1999; Telting et al. 2014) mounted at the 2.56-m Nordic Opti¬ 
cal Telescope (NOT) of Roque de los Muchachos Observatory in 
La Palma, Spain. The observations were performed on 11 August 

2014, 5 November 2014, 20 November 2015, and 26 November 

2015. The data were reduced using standard procedures, which in¬ 
clude bias subtraction, flat fielding, order tracing and extraction, 
and wavelength calibration. The extracted spectra have a S/N ratio 
of 45-55 per pixel at 5500 A. 

Following the same spectral analysis procedure described in 
Rappaport et al. (2015), we use the SPECTRUM code to calcu¬ 
late a grid of synthetic spectra using ATLAS9 models. We then 
use the co-added FIES spectrum to determine the stellar effective 
temperature Tefi, surface gravity log g, projected rotational velocity 
V sini, metal abundance [M/H], and spectral type of KIC 8462852 
(Table 3). The plots in Eigure 5 show selected regions of the ob¬ 
served spectrum (black) along with the best fit model (red). The 
temperature we derive (Teff = 6750 ± 140 K) is consistent with 
the photometric estimate of Teff = 6584^279 K from the revised 
Kepler Input Catalog properties (Huber et al. 2014), as well as 
with Ted = 6780 K derived from the empirical (V — K) color- 
temperature relation from Boyajian et al. (2013). The projected ro¬ 
tational velocity we measure v sin i = 84 ± 4 km s“^ is also well 
in line with the one predicted from rotation in Section 2.1, if the 
0.88 d signal is in fact the rotation period. Overall, the star’s spec¬ 
trum is unremarkable, as it looks like an ordinary early F-star with 
no signs of any emission lines or P-Cygni profiles. Finally, we use 
the stellar properties derived from our spectroscopic analysis to es¬ 
timate a stellar mass M = 1.43 Mq, luminosity L = 4.68 L©, 
and radius R = 1.58 Rq, corresponding to a main-sequence F3 V 
star based on the empirical calibrations from Pecaut & Mamajek 
(2013)^. Combining the radius (assuming a conservative value of 
20 % for the radius error), projected rotational velocity, and rotation 
period (Section 2.1 ), we determine a stellar rotation axis inclination 
of 68 ± 29 degrees. 

While interstellar medium features are not typically related 
to indicators of astrophysically interesting happenings in stars, we 
note the presence of stellar and interstellar Na D lines in our spec¬ 
tra. In the bottom panel of Figure 5, we show a close up of the 
region containing the Na D lines (AA5890, 5896A). Within the two 
broad stellar features, there are two very deep and narrow Na D 
lines with split line profiles, indicating the presence of two discrete 
ISM clouds with different velocities of ~ 20 km s“^. 

2.3 Imaging 

Figure 6 shows the UKIRT image of KIC 8462852 as well as a simi¬ 
larly bright source ~ 40" away. The PSF of KIC 8462852 is asym- 

^ http : //www .pas.ro Chester.edu/-emamajek/EEM_ 
dwarf_UBVIJHK_colors_Teff . txt 



Figure 5. NOT spectrum closeups for KIC 8462852, the best fit stellar 
model shown in red. Panels show region near H «, H /3, Mg, and Na D (top 
to bottom). The bottom panel shows both the stellar (broad) and interstellar 
(narrow) counterparts of the Na D lines. Refer to Section 2.2 for details. 


metric by comparison, leading us to speculate that KIC 8462852 
has a faint companion star about 1.5 — 2" away. 

We observed KIC 8462852 on UT 16 Oct 2014 using the nat¬ 
ural guide star adaptive optics (AO) system (Wizinowich et al. 
2000) of the 10-meter Keck II Telescope on Mauna Kea, Hawaii. 
We used the facility IR camera NIRC2 and the J (1.25 p.m), H 
(1.64 pm), and K (2.20 pm) filters from the Mauna Kea Observato¬ 
ries (MKO) filter consortium (Simons & Tokunaga 2002; Tokunaga 
et al. 2002). We used NIRC2’s narrow camera, which produces a 
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0.00994" pixel”'^ scale and a 10.2" field of view. Conditions were 
cloudy with variable seeing, around 1" FWHM. KIC 8462852 was 
observed over an airmass range of 1.26-1.28. 

The AO-corrected images have full widths at half maxima 
(FWHMs) of 39 mas, 43 mas, and 51 mas at JHK, respectively, 
with RMS variations of about 1-3%. We obtained a series of nine 
images in each filter. The total on-source integration time was 
65 seconds per filter. The images were reduced in a standard fash¬ 
ion using custom scripts written in the Interactive Data Language 
(IDL). We constructed flat fields from the differences of images 
of the telescope dome interior with and without lamp illumination. 
We subtracted an average bias from the images and divided by the 
flat-field. Then we created a master sky frame from the median 
average of the bias-subtracted, flat-fielded images and subtracted it 
from the individual reduced images. The individual reduced images 
were registered and stacked to form a final mosaic (Figure 7). 

As suspected from the asymmetric UKIRT image, the Keck 
AO images reveal an obvious faint companion at a separation of 
1.95" and position angle of 96.6°. To measure the flux ratios 
and relative positions of the two components, we used an analytic 
model of the point spread function (PSF) as the sum of two ellip¬ 
tical Gaussian components, a narrow component for the PSF core 
and a broad component for the PSF halo, as we have done for other 
binaries (Liu et al. 2008). For the individual images obtained with 
each filter, we fitted for the flux ratio, separation, and position angle 
of the binary. To correct for optical distortions in NIRC2, we used 
the calibration of Yelda et al. (2010). The averages of the results 
were adopted as the final measurements and the standard deviations 
as the errors (Table 3). 

It is unclear whether this is a physical or visual binary, though 
given the delta magnitude and separation, the chance alignment of 
the companion being a background or foreground object is only 
~ 1% (Rappaport et al. 2014). At ~ 2% of the flux of the brighter 
star, this would be a ~ 0.4 Mq M2V star, if it is indeed at the 
same distance as our target F star (Kraus & Hillenbrand 2007). The 
JHK colors are also consistent with the companion being a dwarf, 
not a giant (Bessell & Brett 1988). If we take the magnitude of 
KIC 8462852 as V = 11.705, and the absolute visual magnitude 
of an F3V star to be 1/ = 3.08 (Pecaut & Mamajek 2013), then 
the (reddened) distance modulus is 8.625. We derive a de-reddened 
distance of ~ 454 pc using E(B — V) = 0.11 (Section 2.4; cor¬ 
responding to a (/-band extinction of Ay = 0.341). Assuming 
the fainter star is associated with the main F-star target, the angu¬ 
lar separation of ~ 1.95" translates to a distance of ~ 885 AU. 
At this separation, the second star cannot currently be physically 
affecting the behavior of the Kepler target star, though could be af¬ 
fecting bodies in orbit around it via long term perturbations (see 
Kaib et al. 2013). If such a star is unbound from KIC 8462852, but 
traveling through the system perpendicular to our line of sight, it 
would take only 400 years to double its separation if traveling at 
10 km sec~^. So, the passage would be relatively short-lived in 
astronomical terms. 

We also obtained Speckle observations of KIC 8462852 on the 
night of UT 22 Oct 2015 using the DSSI instrument on the WIYN 
telescope located on Kitt Peak (Howell et al. 2011). Observations 
were made simultaneously in two filters with central wavelengths 
at 692 and 880 nm. Both filters show the source to be single, with 
no visible companion observed to within 0.08 arcsec and brighter 
than a delta magnitude of 3.8 and 4.2 magnitudes (for the 692 and 
880 nm filters, respectively). The companion star seen in the Keck 
NIRC2 image was not detected, favoring the conclusion that it is 
an M-dwarf, which would be too faint to be detected in the reddest 



Figure 6. UKIRT image for KIC 8462852 and another bright star for com¬ 
parison, showing that it has a distinct protrusion to the left (east). For ref¬ 
erence, the grid lines in the image are 10" X 10". The color coding is 
logarithmically scaled. Refer to Section 2.3 for details. 



Figure 7. Keck AO iT-band image for KIC 8462852 showing the com¬ 
panion was detected with a 2" separation and a magnitude difference 
A.H = 3.8. The color coding is logarithmically scaled. Refer to Section 2.3 
for details. 

DSSI filter (880 nm). However, it is important to note that these 
speckle results provide an independent confirmation of the results 
from Keck AO: KIC 8462852 has no additional companion down to 
a separation of ~ 20 AU detectable within the relative brightness 
limits with each instrument. 

Finally, we speculate that the 10-20 day periodicity discussed 
in Sect. 2.1 might actually arise on the neighboring faint M star. The 
amplitude of those modulations are ~500 ppm of the total target 
flux. If they arise on the M-star, then their fractional modulation of 
that star would be as high as 3%, which would not be unusual for 
an M star. 

2.4 Spectral energy distribution 

The spectral energy distribution (SED) of KIC 8462852 including 
optical, 2MASS (Skrutskie et al. 2006), (ALL)WISE (Wright et al. 
2010), and Galex NUV (Morrissey et al. 2007) flux densities is 
shown in Figure 8. Optical photometry in BV{RI)c filters was 
obtained by the 90 cm Schmidt telescope of the Konkoly Obser¬ 
vatory at Piszkesteto Mountain Station. For standard magnitudes 
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Wavelength [[im] 

Figure 8. SED for KIC 8462852. The Black solid line is a model for a star 
with Tgfj = 6750 K and E(B — V) = 0.11. Flux calibrated photometry 
are plotted in red, where the extent of the “error-bars” in the X-direction 
indicate the wavelength range of each bandpass and the Y-direction shows 
the error of the flux measurement. Refer to Section 2.4 for details. 

GD391 ABCE photometric standard stars were used as comparison 
(Landolt 2013). Photometric magnitudes are listed in Table 3. 

In order to study whether the system exhibits excess at mid- 
infrared wavelengths, we first fitted an ATLAS9 atmosphere model 
(Castelli & Kurucz 2004) to the photometric data points between 
0.15 and 3.6p.m. From the grid of model atmospheres we selected 
the one that has the closest metallicity, surface gravity, and effective 
temperature to those derived from our spectroscopic study. Thus 
we fixed TeS, log g, and [Fe/H] parameters to 6750 K, 4.0, and 0.0, 
respectively, and only the amplitude of the model and the redden¬ 
ing were fitted. The best fitted photospheric model is displayed in 
Figures. We derive a reddening of 0.11 ± 0.03 mags. By com¬ 
paring the measured W2 and W3 WISF flux densities at 4.6 and 
11.6 pm (at 22 pm we have only an upper limit) with the predicted 
fluxes derived from the photosphere model we found them to be 
consistent, i.e., no excess emission can be detected at mid-infrared 
wavelengths. This lack of significant IR excess is independently 
confirmed using warm Spitzer/IRAC data by Marengo et al. (2015). 

However, this does not exclude the existence of a colder de¬ 
bris disk or a warmer, but relatively tenuous disk. Assuming that 
the emitting grains act like a blackbody, we can derive their char¬ 
acteristic temperature at a specific stellar-centric distance. Using 
this approach, we compute the SFD of a narrow dust belt located 
at a distance of 1, 2, 3, 5, and 10 AU from a star with a luminos¬ 
ity of 4.7 Lq, corresponding to the main-sequence stage (Pecaut 
& Mamajek 2013). The W3 and 1U4 band photometry were then 
used as upper limits to set the amplitude of the excess. Figure 8 
shows the result of these computations and summarizes the fun¬ 
damental disk properties (dust temperature, upper limits for frac¬ 
tional luminosity) of the dust belts at different radii. It is worth 
noting that this very simple model accounts only for large black- 
body grains, smaller (pm-sized) grains are ineffective emitters and 
may be heated to higher temperatures compared to larger grains at 
the same location. We revisit this analysis in more detail later in 
Section 4 (also see Figure 1 1). 

2.5 Ground-based photometric surveys 

We reviewed the ~ 700 photometric intensities from the years 
1900 - 2000 from the Digital Access to a Sky Century Harvard 


(DASCH) project'^ (Grindlay et al. 2012). The error bars on the 
photometry are about ~ 10%. At this level, we found the star did 
not do anything spectacular over the past 100 years. However, if it 
underwent several ~ 20 % dips in flux lasting for several days each 
during that period, the chances are high that there were no plates 
exposed at those times. 

SuperWASP data (Butters et al. 2010) are unremarkable for 
KIC 8462852. We note that there is a 0.2 magnitude offset between 
the available SuperWASP data sets. However, we see the same off¬ 
set when comparing its photometry with a similarly bright source 
nearby KIC 8462852. Thus, we reject this being real (e.g., due to a 
flaring event, etc.). 

Unfortunately, KIC 8462852 falls outside the area covered by 
the KELT network (T. Beatty, private communication). 

2.6 Limits on a close companion 

We use the four FIFS spectra (Section 2.2) to measure the pres¬ 
ence of any Doppler shifts induced by a companion. We traced 
the radial velocity (RV) drift of the instrument by taking long- 
exposed ThAr spectra in a bracketed sequence, i.e., right before 
and after each target observation. RV measurements were derived 
by cross-correlating the target spectra with the rotationally broad¬ 
ened best fitting Kurucz model. The RV measurements are listed in 
Table 2 along with the error bars and the barycentric Julian dates 
in barycentric dynamical time. To within the ~ 400 m s”'^ uncer¬ 
tainties in the RV measurements, the four values we measure are 
quite consistent with no change at all over the 470 day observation 
interval. 

In order to quantify what limits we can set on the mass of 
an hypothetical close companion star, we carried out the following 
analysis. We assumed a circular orbit because there are insufficient 
data points to fit for the parameters in an eccentric orbit. Then, 
for each in a sequence of 4 x 10® trial orbital periods, P, in the 
range of 0.5 to 3000 days, we fit the four RV points with a sine and 
cosine term to represent the orbit and a systemic 7 velocity. From 
this fit we computed the velocity semi-amplitude K and added its 
2-(j uncertainty to establish a conservative upper limit to K. We 
then used the upper limit on K to compute the corresponding upper 
limit on the mass function. Finally, we solved for the upper limit on 
the mass of the hypothetical close companion by taking the mass 
of the F star to be 1.4 Mq, and assuming three different orbital 
inclination angles (30°, 60°, and 90°). The results are shown in 
Figure 9. The spikes are at values of P where the epochs of the four 
RV measurements are commensurate with being at orbital phase 0, 
and the mass constraints are weaker at these periods. For longer 
periods, the density of these spikes diminishes greatly and the lower 
locus of points can be taken as a likely upper limit on the mass of 
any companion. Therefore, we conclude that for periods between 
~30 and 300 days, the mass of any companion is very unlikely to 
exceed that of a brown dwarf. 

Another diagnostic to constrain the nature of the companion 
uses the FT in Figure 2, which shows no sharp, narrow peaks with¬ 
out harmonics (Section 2.1). With this information, a very basic 
limit can be set on a companion from the lack of observed ellip¬ 
soidal light variations (FLVs). The FLV amplitude Aelv is ex¬ 
pressed as: 

Aelv — l-5(Mc/M*)(i7*/a)® sin^ i (1) 

^ http://dasch.rc.fas.harvard.edu/index.php 
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Table 2. TIES RVs of KIC 8462852 


BID 

(- 2450000) 

RV 

[km/s] 

<^RV 

[km/s] 

6881.5183 

4.160 

0.405 

6966.3635 

4.165 

0.446 

7347.3856 

3.817 

0.406 

7353.3511 

4.630 

0.436 


See Section 2.6 for details. 



1 10 100 1000 
Orbital Period (days) 


Figure 9. Upper limits (2-cr confidence) to a hypothetical companion mass 
from the lack of ELVs (red curves) and lack of RV variations on four oc¬ 
casions (blue curves). Each type of constraint is shown for three different 
assumed orbital inclination angles (30°, 60°, 90°); these are marked di¬ 
rectly on the ELV constraint curves, and can be infeiTed from the dashed, 
dotted, and solid curves, respectively, for the RV constraints. These results 
indicate that there are no objects heavier than a super-Jupiter in close orbits 
with Porb $ 2 days, and likely no heavier in mass than a brown dwarf for 
Porb $ 300 days. Refer to Section 2.6 for details. 

(e.g., Kopal 1959; Carter et al. 2011) where M* and P* are the 
mass and radius of the primary, a and i are the semimajor axis and 
orbital inclination, and Me is the mass of a putative companion. 
Rearranging to express a as the orbital period P using Kepler’s 
third law, this equation simplifies to: 

Ablv 3.3 X 10“®(Mc/Mj)(ld/P)^ sin^ i (2) 

where now the companion mass Me is expressed in Jupiter masses 
Mj and P is in days. If ELVs were present, we would have seen 
a peak > 50 ppm for all periods shorter than 4 days (> 0.25 cy¬ 
cles day“^) in the FT (Figure 2). 

The limits on the companion mass that we can set from the 
lack of ELVs, as a function of orbital period, are illustrated in Fig¬ 
ure 9. They are plotted as red lines for three different assumed incli¬ 
nation angles. Note that an angle of 90° is not allowed or we would 
have seen (regular) transits; it is shown in this figure for instruc¬ 
tive purposes only. These ELV mass constraints are superposed on 
those discussed above based on the lack of RV differences among 
our four measurements (blue curves). Taken together, these results 
indicate that there is not likely to be a close companion to the F star 
more massive than a super-Jupiter with Porh tS 2 days, nor more 
massive than a brown dwarf for Porb < 300 days. 


Table 3. Properties of KIC 8462852 


Property 

Value 

Method/Reference 

RA (deg) 

301.564392 

KIC 

DEC (deg) 

44.456875 

KIC 

Kp (mag) 

11.912 

KIC 

B (mag) 

12.262 4: 0.008 

90 cm Schmidt (§ 2.4) 

V (mag) 

11.705 4:0.017 

90 cm Schmidt (§ 2.4) 

Be (mag) 

11.356 4:0.024 

90 cm Schmidt (§ 2.4) 

Ic (mag) 

11.051 4: 0.098 

90 cm Schmidt (§ 2.4) 

J (mag) 

10.763 4:0.021 

2MASS 

H (mag) 

10.551 4: 0.019 

2MASS 

K (mag) 

10.499 4:0.020 

2MASS 

W1 (mag) 

10.425 4:0.023 

(ALL)WISE 

W2 (mag) 

10.436 4:0.020 

(ALL)WISE 

W3 (mag) 

10.591 4: 0.123 

(ALL)WISE 

W4 (mag) 

9.423“ 

(ALL)WISE 

Rotational period (d) 

0.8797 4: 0.0001 

FT(§ 2.1) 

Spectral type 

F3 V 

Spectroscopy (§ 2.2) 

Teff (K) 

6750 4: 120 

Spectroscopy (§ 2.2) 

logs (ogs) 

4.0 4: 0.2 

Spectroscopy (§ 2.2) 

[M/H] (dex) 

0.00 4: 0.10 

Spectroscopy (§ 2.2) 

V sin i (km s“^) 

84 4:4 

Spectroscopy (§ 2.2) 

distance (pc) 

454 

Distance modulus (§ 2.3) 

E{B - V) (mag) 

0.11 4: 0.03 

SED (§ 2.4) 

Binary separation (arcsec) 

1.96 

Keck AO (§ 2.3) 

Binary position angle (deg) 

96.6 

Keck AO (§ 2.3) 

A J (mag) 

4.209 4: 0.044 

Keck AO (§ 2.3) 

AH (mag) 

3.840 4: 0.017 

Keck AO (§ 2.3) 

AK (mag) 

3.619 4:0.012 

Keck AO (§ 2.3) 


“Upper limit. 


2.7 Space motion and age 

Using our distance estimate of 454 pc (Section 2.3), the radial ve¬ 
locity obtained from the FIES spectrum (Section 2.6), and proper 
motions and positions from the UCAC4 catalogue we computed 
the Galactic space motion of the target, yielding 4-31.5, —2.5, and 
4-10.2 km s”'^ for the U (defined as positive toward the Galac¬ 
tic center), V, and W velocity components, respectively. Young 
disk population stars have low velocity dispersion and they occupy 
a special region within the velocity space. Based on the studies 
of Eggen (1989), Leggett (1992) defined a box by —50 < U < 
-420 km s-\ -30 < V < 0 km s"\ and -25 < W < 10 km s"\ 
which includes most of the young disk stars in our neighborhood. 
Our target lies outside of this box. In fact, its galactic space mo¬ 
tion - especially the U component - deviates significantly from 
the characteristic space motion of any nearby young (< 100 Myr) 
kinematic groups, open clusters, and star forming regions (Makarov 
2007; Mamajek 2015). Altogether, it implies that KIC 8462852 
likely does not belong to the youngest stellar population. 

In making this distance estimate, we assumed that 
KIC 8462852 is a main-sequence star (Section 2.3). We note 
that assuming a pre-main or post-main sequence phase does not 
change our previous conclusion. These evolutionary stages would 
be accompanied by larger luminosities and thereby larger dis¬ 
tances. This would result in a galactic space motion that deviates 
even more significantly from that of typical young disk stars. 
Unfortunately, our star falls outside the region where empirically 
calibrated age diagnostics such as chromospheric activity or stellar 
rotation period can be used (e.g., Mamajek & Hillenbrand 2008). 
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2.8 Similar dippers in the Kepler field? 

The anomalous dips in KIC 8462852 were serendipitously found 
by the Planet Hunter citizen science group. Due to its aperiodic 
nature, it likely never would have been flagged/recovered by most 
searches for transits, eclipsing binaries, or asteroseismologically in¬ 
teresting stars. However, knowing the existence of KIC 8462852’s 
light curve, we naturally wondered if there are, in fact, numerous 
other such objects in the main-field Kepler data base. We therefore 
applied a simple algorithm to search the data base for other systems 
similar to KIC 8462852. The algorithm consisted of searching for 
dips with depths of greater than 10% (i.e., normalized fluxes of 
< 0.9) that consist of 5 or more consecutive Kepler long-cadence 
samples (i.e, lasting more than ~ 2.5 hours). In all, this search 
turned up more than a thousand targets with this signature. The 
vast majority of them, however, were due to (1) eclipsing binaries, 
(2) the rotation signature of large amplitude starspots, and (3) some 
obvious Kepler data artifacts. We carefully examined the remain¬ 
ing small number of systems by eye, but could identify none that 
was reminiscent of KIC 8462852. We also lowered the threshold for 
dips to 5%, but the search likewise turned up no candidates that one 
would believe closely resemble KIC 8462852. Of course, some of 
the visual comparison work is necessarily qualitative, but we were 
satisfied that there are at most a few similar systems to be found in 
the main Kepler field. 


3 POSSIBLE EXPLANATIONS OE THE OBSERVED 
DIPPING EVENTS OBSERVED IN KIC 8462852 

The main issue in explaining the peculiar light curve for 
KIC 8462852 is related to the presence of multiple dimming events, 
that are not periodic and of which the D800 single event has a 
smooth, yet highly asymmetric, profile, and the D1500 events are 
the deepest and most complex. Here, we introduce several scenar¬ 
ios to explain KIC 8462852 and discuss how the observational data 
do and do not support each theory. 


3.1 Instrumental effects or data reduction artifacts? 

The Kepler light curve for KIC 8462852 is unique, and we have 
thoroughly explored the raw data for defects/instrumental effects, 
which could cause the observed variations in KIC 8462852’s flux. 
We use the PyKe software tools for Kepler data analysis to check 
the data for instrumental effects. We check the following possibili¬ 
ties: 

• We checked that the same flux variations, i.e., the ‘dips’, are 
present in the SAPJFLUX data set. 

• We verified that data gaps and cosmic rays events® do not co¬ 
incide with the dipping events, as they are prone to produce glitches 
in the corrected fluxes. 

• We verified at the pixel-level that there are no signs of peculiar 
photometric masks used in making the light curves. 

• We verified at the pixel level that the image light centroid does 
not shift during the ‘dipping’ events 

• We inspected the Kepler light curves of neighboring sources 
and find that they do not show similar variability patterns in their 
light curves. 


® The times of these events are recorded in the headers of the fits files 


• We determined that CCD cross talk, reflection, and column 
anomaly cannot be the cause (Coughlin et al. 2014). 

• We verified with the Kepler team mission scientists that the 
data were of good quality. 

This analysis concludes that instrumental effects or artifacts in 
the data reduction are not the cause of the observed dipping events, 
and thus the nature of KIC 8462852’s light curve is astrophysical 
in origin. 

3.2 Intrinsic variability? 

An example of a class of stars which display intrinsic variability 
are the R Coronae Borealis (RCB) type variables. These are highly 
evolved F-G supergiants (e.g., Clayton 1996) that have light curves 
which show pulsations (on the order of months) and irregular deep 
dips (lasting weeks to months). Their “dipping” variability is asso¬ 
ciated with formation of clouds that obscure the photosphere, and is 
often observed as a sharp decrease in flux followed by a more grad¬ 
ual, and sometimes staggered, recovery. In the case of KIC 8462852 
the time scales of the dips are different than those of a RCB vari¬ 
able. Likewise, the ingress at D800 has a gradual decrease in flux, 
which is inverse to what is expected in a RCB, and the dip shapes 
at D1500 are also non-characteristic of a RCB. Lastly, the spectro¬ 
scopic log g and v sin i are far from those of a supergiant. These 
items together strongly rule out the possibility of KIC 8462852 be¬ 
ing a RCB variable. 

Another possibility is the self-emission of disk material from 
the star itself, as in the case of Be-stars. Be stars are rapidly rotating 
(almost near breakup) stars that are usually of spectral class O and 
B, but sometimes A, and exhibit irregular episodic outbursts. Usu¬ 
ally these outbursts are in emission, but in some cases it can also 
result in dimming (see Hubert & Floquet 1998). Be stars also often 
exhibit quasi-periodic oscillations in the range of ~ 0.5 — 1.5 days. 
This also fits the bill for what we see in the FT of KIC 8462852 
(§ 2.1). It has been hypothesized (e.g., Rappaport & van den Heuvel 
1982) that most, if not all, Be stars have a binary companion which 
originally transferred mass to the current Be star to spin it up to 
near breakup (the remnant of that star is sometimes found to be a 
neutron star). The periods of these binaries range from a couple of 
weeks to thousands of days (perhaps longer). If KIC 8462852 is a 
Be star, we would get an unprecedented look into the inner disk be¬ 
havior. In such as case, the broad peak in the FT at frequencies just 
below the 0.88 d periodicity could be explained by ejected mate¬ 
rial in a so-called “excretion disk” that is moving outward but with 
roughly Keplerian velocity. 

The lack of observed IR excess does not support the existence 
of an excretion disk. There is also an absence of Ha emission in 
the star’s spectrum (although, as noted above. Be star Ha emission 
is known to be variable and turn off and on with timescales from 
days to years). Furthermore, the temperature of KIC 8462852, Tes 
= 6750 K, is too cool to be a Be star. It is also unlikely to have been 
spun-up by a donor star whose remnant is still orbiting the F star 
because of the constraints set by the four RV measurements and the 
limits on any ELVs (see Section 2.6). Though, we cannot rule out 
remnants orbiting with P > a few years. 

3.3 Extrinsic variability? 

3.3.1 Related to the secondary star 

We first consider whether KIC 8462852’s flux is contaminated by 
the nearby M-dwarf detected with high-resolution images (§ 2.3). 
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Whether or not the system is hound, the faint companion con¬ 
tributes light in the Kepler photometric aperture, which in turn af¬ 
fects the observed signal in the light curve. Our observations show 
that the flux ratio in the infrared is ~ 30, which translates to a 
factor of several hundred in the Kepler bandpass. Thus, the max¬ 
imum imprint that the M-dwarf has on the light curve variabil¬ 
ity is ~30 mmags; this is insufficient to make an impression on 
KIC 8462852’s light curve at anything greater than ~3%, and, in 
particular, it couldn’t possibly explain any of the large dips. 

3.3.2 Occultation by circumstellar dust clumps 

The dips could be readily explained in terms of occultation by an 
inhomogeneous circumstellar dust distribution. However, this does 
not mean that the dust distribution that would be required to explain 
the observations is physically plausible or would necessarily apply 
to KIC 8462852. 

Inhomogeneous dust distributions have been invoked to ex¬ 
plain dips seen towards some young stars such as UX Orionis 
or AA Tau-like “dipper” systems (Herbst et al. 1994; Herbst & 
Shevchenko 1999; Morales et al. 2009; Cody et al. 2014; Ansdell 
et al. 2015). At an age of only a few tens of Myrs, these dipper 
stars have V -band light curves characterized by sporadic photomet¬ 
ric minima with amplitudes of 2 - 3 magnitudes and with durations 
of days to many weeks. These objects also generally exhibit strong 
infrared excess, starting at ~ 2 — 5 pm and show signs of accretion 
(emission) in their spectra. However, in contrast to such systems, 
KIC 8462852 has no detectable IR excess or accretion signature to 
suggest that it is a young T Tauri star (Sections 2.2, 2.4). Thus a sce¬ 
nario in which material in a gas-dominated protoplanetary disk oc¬ 
cults the star due either to accretion columns or non-axisymmetric 
azimuthal or vertical structure in the inner disk (e.g. Herbst et al. 
1994; Herbst & Shevchenko 1999; Bouvier et al. 1999; McGinnis 
et al. 2015) is strongly disfavoured. 

We therefore are left to consider scenarios that could arise 
around a main-sequence or weak-line T Tauri star that has dispersed 
its protoplanetary disk, but still hosts a gas-poor planetary system 
that may include planets, asteroids, and comets. The “clumps” of 
dust passing in front of the star could perhaps lie within an opti¬ 
cally thin asteroid belt analogue that is otherwise undetected, or 
be more isolated objects such as remnants of a broken up comet. 
As in the above scenarios, the typical minimum sizes of the dust 
grains are ~pm. (e.g. Backman & Paresce 1993), which are able to 
cause stellar variation by absorbing and scattering starlight at opti¬ 
cal wavelengths. Before considering such scenarios in more detail, 
we start with some scenario-independent constraints that can be 
gleaned from the observations. 


4 SCENARIO-INDEPENDENT CONSTRAINTS 

To understand what could be the origin of the clumps it would help 
to know where they are located in the system, how big they are, and 
how long they last. To aid with this discussion. Figure 10 shows 
some scenario-independent constraints on the size and orbital dis¬ 
tance of the clumps that are discussed further below. The only as¬ 
sumption for now is that the clumps are on circular orbits, but this 
assumption is relaxed later in Section 5.4. Some of the constraints 
also assume the clumps to be opaque, but again this assumption is 
relaxed later. 

Dip duration: The timescale fdip for the transit of a clump of 
radius s with transverse velocity Vt across the equator of a star with 



0.1 1.0 10.0 100.0 
Semimajor axis (au) 


Figure 10. Size vs. semi-major axis parameter space for optically thick, 
spherical dust clumps on circular orbits around a star of M, = 1.43 Mq 
and R* = 1.58 Rq. The solid lines represent dips of equal duration (as 
labelled). Dotted lines show minimum clump sizes for dips of different 
depths. Vertical dashed lines show where the orbital period is 1500 days, 
and where the light curve gradient for an optically thick “knife edge” could 
be as high as 0.5 d“^. Diagonal dashed lines show Hill radii of planetes- 
imals of different sizes, assuming a density of 3 g cm“®. Combined, the 
period, gradient, and duration constraints in the circular orbit scenario sug¬ 
gests the clumps lie between ~ 3 — 10 AU, and have sizes similar to the 
star. 


radius i?* is fdip = 2 (s + R*) jvt. If the clump is on a circular 
orbit around a star of mass M* with semi-major axis a, and is much 
less massive than the star, then 

s« 1.85fdip (j -R*, (3) 

for a is in units of AU, M* in Mq, s and i?* in Rq, fdtp in days. 
Thus, the several-day duration of the events for KIC 8462852 sug¬ 
gests that the clumps are either close-in and large compared to the 
star, or far-away from the star and small. However, clumps that are 
too distant move too slowly across the stellar disk to explain the ob¬ 
served duration regardless of their size; e.g., a 3-day duration dip 
cannot arise from a clump beyond ~ 15 AU. 

Dip depth: A minimum clump size is set by the depth of the 
dimming events, which we characterise as 1 minus the normalised 
flux, which we call r. For example, even if the clump is completely 
opaque, the maximum dip depth is max(T) = (s/R*)^. The deep¬ 
est r = 20% dimming event at D1500 thus implies that at least 
some clumps are a sizeable fraction of the stellar size. A dip caused 
by a fully optically thick symmetrical clump would also have a 
characteristic symmetrical shape which does not resemble those 
observed (e.g., panel ‘c’ in Figure 1), so this can be regarded as 
a strong lower limit. While there appear to be a range of event du¬ 
rations, the duration of the deepest events is at most about 3 days. 
The middle solid line in Figure 10 (for fdip = 3 d and a depth of 
r = 20%) therefore decreases the outer limit on the clump loca¬ 
tions mentioned above to closer to 8 AU. 

Light-cun’e gradient: A similar, but independently derived, 
outer constraint on the clump location can be set by examining the 
gradients in the light-curve, which are at most half of the total stel¬ 
lar flux per day (i.e. 0.5 d“^ when the light curve is normalised 
to 1). Orbiting material can change the light-curve most rapidly 
when it is optically thick and passing the stellar equator (i.e., the 
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Figure 11. Fractional luminosity limits (blue lines) and an estimate of the 
system dust content from the light curve (green line). The dust level is con¬ 
strained to lie below the blue line by the WISE photometry (4.6 /rm, 12 fim, 
and 22 fim). The green line integrates the optical depth in the light curve 
assuming that clumps are similar in size to the star and on circular orbits. If 
the clumps lie beyond about 0.2 AU the IR non-detection of the dust is un¬ 
surprising, although many scenarios require more emission than that from 
dust seen to pass along our line-of-sight to the star. Refer to Section 4 for 
details. 


“knife edge” model of van Werkhoven et al. 2014). The high rate of 
change in the KIC 8462852 light curve translates to a lower limit on 
the transverse velocity of the orbiting material of about 9 km s“^, 
which corresponds to an upper limit of 13 AU for material on cir¬ 
cular orbits, although as discussed in Section 5.4, this upper limit 
is closer to the star if the clump is optically thin. 

Non-periodicity: The lack of evidence for a clear periodicity in 
the dips in the observed light-curve excludes orbital periods shorter 
than ~ 1500 days, which thus constrains the location to he beyond 
about 3 AU. This constraint could be broken if the clumps disperse 
within a single orbit. Likewise, if the two deep dipping events at 
D800 and D1500 are from the same orbiting body (or bodies), a 
period of 700 - 800 days remains a possibility. 

Gravitational binding: To address the survival of the clumps, 
we note that in any scenario where the clumps are not self- 
gravitating, they cannot be long-lived in the face of orbital shear 
(e.g. Kenyon & Bromley 2005) and their internal velocity disper¬ 
sion (e.g. Jackson & Wyatt 2012). Figure 10 therefore shows plan- 
etesimal sizes required to retain dust clouds within their Hill sphere, 
F?Hiii = a(Mpi/[3M,])^^®, as one way of ensuring long-lived 
clumps. 

Thus, under the assumption of circular orbits, the depth, dura¬ 
tion and lack of periodicity of the dimming events constrains their 
location to a region roughly corresponding to that occupied by the 
giant planets in the Solar System (i.e., between the green dashed 
lines). Clump sizes would thus be comparable to, but larger than, 
the star (i.e., above the uppermost horizontal dotted purple line), 
and they would have to have high, but not necessarily unity opti¬ 
cal depth. It might be possible to explain the clumps as dust bound 
to planetesimals larger than around 1000 km, which means such 
planetesimals are not necessarily large enough for direct transit de¬ 
tection (the lack of which could provide another constraint). 

Infrared excess: Another constraint on the origin of the 
clumps comes from the lack of infrared emission (Section 2.4). 



Time (days) 


Figure 12. Inverted light curve for KIC 8462852 portraying the blocking 
factors needed to reproduce the light curve as a function of time. Refer to 
Section 4 for details. 


Assuming the clumps are larger than the star, the Kepler light 
curve provides blocking factors needed as a function of time. 
In (normalized flux), where ln(normalised flux) « r for small 
r, as shown in Figure 12. This optical depth and the assumption 
that the clump crosses the star at its orbital velocity allows conver¬ 
sion to optical depth as a function of distance along the clump. The 
dimming events therefore allow an estimate of the minimum pos¬ 
sible cross-sectional area atot of dust in orbit around the star. That 
is, 

cTtot = Vth J r{t)dt, (4) 

where the light-curve yields f T{t)dt ~ 0.86 days, Vt is the ve¬ 
locity of the clumps (assumed to be uniform at circular velocity 
for a given semi-major axis), and h the “height” of the clumps (i.e. 
their size along the dimension perpendicular to their velocity). The 
height of the clumps is assumed to be 2 R,, though it could be 
higher if not all of the clump crosses the stellar disk (e.g., this could 
be assumed to be 7rs/2 for large spherical clumps passing directly 
across the star). This calculation gives the minimum possible cross- 
sectional area as 

atot = 2.6 X (5) 

where a is in AU, the dependence on which arises from the velocity 
at which the clump crosses the star. 

This cross-sectional area can then be converted to frac¬ 
tional luminosity at a given distance from the star using / = 
atot/(4TTa^). The blue lines in Figure 11 show the limits on the 
dust fractional luminosity / = Ldust/i* derived from the SED 
(Section 2.4). These can be thought of as the maximum luminos¬ 
ity of blackbodies at a range of dust temperatures (or stellocentric 
radii) that fit under the WISE photometry. The dust estimate from 
Equation 4 is shown as a green line, and the fact that it lies below 
the blue line at all radii beyond 0.2 AU indicates that it is perhaps 
not particularly surprising that no mid-IR excess was seen. 

However, this dust area estimate is only a lower limit since 
it only includes the dust which passed in front of the star during 
the lifetime of the Kepler mission. The true area would be larger 
if there are more clumps further along the orbit which have yet to 
pass in front of the star, and could also be larger if the dips do not 
capture all of the cross-sectional area in their clumps. Furthermore, 
for some specific scenarios discussed in the following sections, the 
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presence of clumps that pass in front of the star requires the exis¬ 
tence of other clumps that do not pass along our line-of-sight. The 
lack of infrared emission thus places constraints on how many such 
clumps there are in the system. For example, Figure 11 shows that 
for clumps at a few AU the cross-sectional area can only be in¬ 
creased by 3 orders of magnitude before it is detectable by WISE. 
The calculation is further complicated should the clumps be con¬ 
sidered to be short-lived, or on non-circular orbits. 

Mass estimates: The minimum possible cross-sectional area 
required to cause the observed dips, atot (Equation 5), can also 
be used to determine a minimum possible dust mass, rritot- If the 
dust all has the same diameter D and density p then rritot/otot = 
2pD jZ, resulting in a total mass of 6.7 x 10^® g for 1 pm diameter 
dust of density 3 g cm“® orbiting at 3 AU (and scaling as a~^^^pD 
for different assumptions). If all of this mass were put in a single 
body of the same density p, this would have a diameter of 16 km. 
This illustrates that the minimum mass of the parent body required 
to cause this phenomenon is approaching the mass of comet Hale- 
Bopp. However, this calculation has a few caveats. For one, the 
value derived for atot is an absolute minimum given that it only 
accounts for the material which passed in front of the star during 
the observations. It also does not account for the possibility that the 
dust in the clump has a range of sizes. For example, for dust with 
a power law size distribution with index of 3.5 (Dohnanyi 1969) 
extending from Dmin to Dmax, the ratio of mass to cross-sectional 
area scales oc p\/DtninDmox- Thus, this estimate would be 100 
times larger than that derived above if the size distribution extended 
from 1 pm up to 1 cm. 

Given these basic constraints we now consider several scenar¬ 
ios that may explain the observations. The first two are related to 
collisions within an asteroid belt (Section 5. 1) or unstable planetary 
system (Section 5.2), the third considers dust that orbits within the 
Hill spheres of large planetesimals which may reside in an asteroid 
belt but are not required to collide (Section 5.3), and the fourth is 
that the dips are the passage of a series of fragments from a broken- 
up comet or asteroid on a highly elliptical orbit (Section 5.4). 


5 SPECIFIC OCCULTATION SCENARIOS 

5.1 Aftermath of catastrophic collisions in asteroid helt 

One possibility is that the dimming events are caused by dust 
thrown off in collisions between planetesimals in an otherwise un¬ 
seen asteroid belt analogue (e.g., Wyatt & Dent 2002; Zeegers et al. 
2014). The dust clouds created in these destructive collisions ex¬ 
pand at roughly the planetesimals’ escape velocity from the collid¬ 
ing bodies, eventually spreading and shearing out to form a smooth 
dust component in which the clumps reside. Such a scenario is a 
promising explanation for the star RZ Psc (de Wit et al. 2013), 
though in that case evidence that the underlying asteroid belt ex¬ 
ists is given by a strong IR excess. 

There are several problems with this scenario as applied to 
KIC 8462852 however. Probably the most fundamental of these is 
the absence of an IR excess from the smooth component. This is be¬ 
cause for every clump we see, remembering that these were inferred 
to be slightly larger than the star, there should be many more that 
have spread out. The infrared emission from the dispersed clumps 
would likely sum up to a detectable level, even before counting 
dust produced in non-dip forming events. Moreover we should see 
dips from the clumps in the middle of being dispersed (i.e., dips 
with longer duration albeit lower optical depth), as well as dips 


with a continuum of depths and durations from the many different 
scales of planetesimal impacts that would occur. The clustering of 
dips at D1500 also points to these events being correlated which 
is hard to reconcile with this scenario, though the planetesimals in 
the belts could be shepherded by planets into confined azimuthal 
regions (e.g., Wyatt 2003; Nesvomy et al. 2013). 

5.2 Aftermath of giant impact in planetary system 

A possible way around the issues in Section 5.1 is to invoke dust 
thrown off in a single collision, perhaps analogous to the Earth- 
Moon system forming event (Jackson & Wyatt 2012). In this case 
there need not be an underlying asteroid belt, as the collision could 
be between planets whose orbits recently became unstable, or be¬ 
tween growing planetary embryos. Such events are expected to re¬ 
sult in strong IR excesses (e.g. Jackson & Wyatt 2012; Genda et al. 
2015), and are indeed seen in systems such as HD 172555 where 
giant impacts are the favored explanation (Lisse et al. 2009). In this 
scenario, the putative collision would need to have occurred be¬ 
tween the WISE observation taken in Kepler Q5 and the first large 
dip at D800. The dip at D1500 is then interpreted as the same ma¬ 
terial seen one orbit later, with the ~ 750 day period implying an 
orbit at ~ 1.6 AU. The difference in the dip structure from D800 
to D1500 could arise because the clump(s) created in the original 
impact are expanding and shearing out. This scenario therefore pre¬ 
dicts that KIC 8462852 may now have a large mid-IR excess, but 
the most recent IR observations taken in 2015 January with Spitzer 
IRAC show no significant excess for KIC 8462852 (Marengo et al. 
2015). However, non-detection of an excess would not necessarily 
rule this scenario out, as the dust levels derived in Section 4 (which 
account for the dust seen passing in front of the star) were shown 
to be consistent with a non-detection. A more robust prediction is 
that future dimming events should occur roughly every 750 days, 
with one in 2015 April and another in 2017 May. 

Two new issues arise with this scenario however. Eirstly, if the 
period of the orbiting material is a few years, what is the origin of 
the two small 0.5% dips seen in the first few hundred days (D140 
and D260; Table 1), and why did they not repeat 750 days later? 
It is a concern that these could require the existence of an outer 
planetesimal belt, which may contradict the lack of infrared emis¬ 
sion to this star. Perhaps more problematic is the probability that 
this star (of unknown age) should suffer such an event that occurs 
within a few-year window between the WISE observation and the 
end of the prime Kepler mission, and that the geometry of the sys¬ 
tem is such that material orbiting at ~1.6 AU lies almost exactly 
between us and the star. Taking this few year window, the main 
sequence lifetime, and an optimistic estimate for the scale height 
of giant impact debris, and the number of Kepler stars observed, 
this suggests that every star would have to undergo 10'* such im¬ 
pacts throughout its lifetime for us to be likely to witness one in 
the Kepler field. Thus, while this scenario is attractive because it 
is predictive, the periodicity argument may be inconsistent, and the 
probability of witnessing such an event may be very low (though 
of course difficult to estimate). 


5.3 Dust-enshrouded planetesimals 

Scenarios in which the clumps can be long-lived are attractive be¬ 
cause they suffer less from being improbable. Thus, one possibil¬ 
ity is that the clumps are held together because they are in fact 
themselves orbiting within the Hill sphere of large planetesimals. 
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They can therefore be thought of as planetesimals enshrouded by 
near-spherical swarms of irregular satellites, which are themselves 
colliding to produce the observed dust. This scenario is therefore 
analogous to that suggested for the enigmatic exoplanet Fomal- 
haut b (Kalas et al. 2008; Kennedy & Wyatt 201 1), which borrows 
from the irregular satellites seen in the Solar System (e.g. Jewitt & 
Haghighipour 2007; Bottke et al. 2010). This scenario suffers from 
several problems. First, the observed dips already require multiple 
large planetesimals. Unless these all orbit within the same plane to 
a high degree (i.e., to within a few stellar radii), there must be many 
more large planetesimals which never (or have yet to) pass in front 
of the star. Debris disks with low levels of stirring are theoretically 
possible (Fleng & Tremaine 2010; Krivov et al. 2013). Flowever, 
these low stirring levels require the absence of large planetesimals 
which through mutual interactions would stir the relative velocities 
to their escape speeds. This is in addition to the problem of filling 
the Hill sphere of such planetesimals almost completely with dust. 
This may be reasonable if the planetesimals are embedded in a belt 
of debris. However, that would incur the problem of the lack of in¬ 
frared excess. The question also remains why the D1500 events are 
so clustered, and why there are several deep dimming events and 
no intermediate ones. A population of planetesimals should have 
a variety of inclinations with respect to our line of sight, so they 
should pass in front of the star at a range of impact parameters and 
cause a range of dip depths. 

A related scenario is that the planetesimals are surrounded 
by large ring systems, similar to that invoked to explain the ~50 
day dimming event seen for ISWASP J140747.93-394542.6 (nor¬ 
mally called “J1407”, Mamajek et al. 2012; van Werkhoven et al. 
2014; Kenworthy & Mamajek 2015). In that case however, a sin¬ 
gle relatively time-symmetric dimming event was seen, whereas 
KIC 8462852 has multiple asymmetric events. Thus, a single ringed 
planet(esimal) would not reproduce the observed light-curve, and a 
scenario with multiple ringed-planetesimals would be essentially 
the same as the irregular satellite scenario above. 

5.4 A family of objects on a comet-like trajectory 

One of the scenario independent constraints considered in Section 4 
was the presence of light-curve gradients as large as 0.5 d~^, which 
results in an upper limit of 13 AU for the clumps’ semi-major axis 
assuming optically thick clumps (Figure 10). However, the star is 
never completely occulted, so this estimate should be corrected for 
the optical depth of the clump r. That is, the steepness of the gra¬ 
dient is diluted either by flux transmitted through a large optically 
thin clump (or by unocculted parts of the star for an optically thick 
small clump). Assuming r = 0.2 the velocity estimate given by 
the gradients is then 5 times higher than assumed in Section 4; 
this would predict a more realistic minimum transverse velocity 
of ~50 km s“^ to cause the observed gradient, which for a circular 
orbit yields a maximum semimajor axis of a = 0.5 AU. While this 
estimate is uncertain, for example because of the unknown optical 
depth structure of the different clumps, this highlights the possibil¬ 
ity that the material may be moving so fast that the velocity for a 
circular orbit is inconsistent with the non-repetition of the events. 

One solution to this problem is that the orbits need not be cir¬ 
cular. That is, we could be seeing material close to the pericenter of 
a highly eccentric orbit, reminiscent of comets seen in the inner So¬ 
lar System at pericenter (Marsden 1967; Sekanina 1984). Comets 
around other stars have also been detected, the first of these being 
found around Beta Pictoris (Smith & Terrile 1984; Lagrange-Henri 
et al. 1989; Beust et al. 1990). We therefore envision a scenario in 
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Figure 13. Size vs. pericenter parameter space for high eccentricity comet¬ 
like orbits. Dotted lines show lower limits on the clump sizes from the dip 
depths. The dashed line is the outer limit set by the light curve gradient, 
noting that this limit decreases with decreasing optical depth, e.g., the limit 
would be at a pericentre that is 25 times smaller than that plotted if the 
clumps have optical depth of 0.2 (line not shown in figure). The dot-dashed 
line is where the clump radius equals the pericenter distance, though the 
clumps could exist above here if they are elongated along the orbital direc¬ 
tion. The solid lines are of constant dip duration. 


which the dimming events are caused by the passage of a series 
of chunks of a broken-up planetesimal on a comet-like orbit. That 
planetesimal may have been analogous to what we refer to in the 
Solar System as a ‘comet’, in which case it could be volatile-rich 
and may have broken up as a result of thermal processes. However, 
it may alternatively have closer analogy with Solar System aster¬ 
oids in having a more refractory composition, which might require 
non-thermal processes such as tidal disruption to break it up. The 
disruption mechanism and composition of the planetesimal are not 
defined for this scenario, just its orbit which is comet-like, and so 
we refer to it here-on as a “comet-like” without bias to their ori¬ 
gin or physical make-up. Regardless of its disruption process, the 
resulting chunks would have to have since spread around the orbit, 
and may be continuing to fragment to cause the erratic nature of the 
observed dips. 

To assess this scenario. Figure 13 revisits the clump - orbit 
parameter space of Figure 10 (discussed in Section 4), but now 
uses the pericenter of the clump’s orbit instead of its semi-major 
axis. The orbits are assumed to be highly eccentric (e ~ 1), with 
the dips arising from material close to pericenter, so that their or¬ 
bital velocity is roughly v/2 times the circular Keplerian velocity 
at that distance. The limits from the dip depths and light-curve gra¬ 
dient are again shown, as are lines of constant dip duration. The 
planetesimal Hill radius lines are not shown, because they are not 
applicable to the cometary scenario considered here, though these 
would be slightly modified versions of those in Figure 10 (see eq. 
B5 of Pearce & Wyatt 2014). In general, the main change com¬ 
pared with Figure 10 is that the higher orbital velocity relaxes the 
constraints on how far out the clumps can be orbiting. However, 
as mentioned above, if the clumps are optically thin (as opposed 
to optically thick as assumed in Figure 13) the constraint from the 
light curve gradient may be more stringent. For example, decreas¬ 
ing the optical depth to 0.2 would result in a transverse velocity 
of 50 km s“^ (see above), thereby moving the light curve gradient 
constraint on the upper limit from 26 AU closer to 1 AU. 


© 0000 RAS, MNRAS 000, 000-000 








KIC 8462852 - Where’s the flux ? 15 


The proximity of the comet-like clump to the star when caus¬ 
ing the dip does not present a problem for this scenario, as it did 
when the clump was on a circular orbit. This is because the peri- 
center distance does not necessarily bear any relation to the period 
with which the comet-like fragments return to pass in front of the 
star. That period is set by the semimajor axis which has the same 
constraint as shown on Figure 10, and there is no such constraint on 
the pericenter in Figure 13. Thus the point of note from Figure 13 is 
that the pericenter could be significantly within 1 AU. Closer peri- 
centers are favored both because this geometry results in a higher 
probability of the clumps occulting the star along our line-of-sight, 
and because of the greater opportunities for fragmentation of the 
bodies. The temperatures of comets (i.e., with volatiles) at such 
close proximity to the star (> 410 K) would render them suscepti¬ 
ble to thermal stresses. The existence of multiple super-Earth plan¬ 
ets orbiting < 1 AU from many main sequence stars also points 
to the possibility that the body could have been tidally disrupted 
in a close encounter with one such planet. It is even possible that 
the body came close enough to the star for tidal disruption in the 
absence of other considerations; e.g., a comet similar to Halley’s 
comet would fall apart by tidal forces on approach to within 3-7 
stellar radii (0.02 - 0.05 AU). By contrast, a rocky body would re¬ 
quire a closer encounter to tidally disrupt. 

For close pericenters it is important to point out that while the 
constraint is discussed in terms of the clump’s radius, the clump 
can not in fact be spherical at that size. Figure 13 shows a blue dot- 
dashed line where the “clump radius” is the same as the pericenter 
distance. At such proximity, the clump could not be elongated in the 
radial direction, but could only be elongated azimuthally along the 
orbit. In fact, this mostly linear clump structure is the correct way 
to visualize debris from the breakup of a comet or planetesimal. 
The small velocity kicks (from fragmentation or tidal disruption) 
would cause a small dispersion in semimajor axis for material in 
the clump, and the resulting differential orbital motion causes the 
material to spread around the orbit. These small kicks do not sig¬ 
nificantly change the periastron distance or the orbital inclination 
angle. 

This scenario is attractive, because comets are known in the 
Solar System to have highly eccentric orbits and disrupt for various 
reasons near pericenter, and infalling comets are the most robust 
explanation for the falling evaporating body (FEB) phenomenon 
seen around many nearby A-type stars (e.g. Kondo & Bruhweiler 
1985; Beust et al. 1990; Welsh & Montgomery 2013; Kiefer et al. 
2014). Also, since fragments of the comet family would all have 
very similar orbits, this mitigates the problem noted in Section 5.1 
that the detection of multiple transits may require orders of mag¬ 
nitude more clumps to be present in the system. Instead, the ob¬ 
served clumps may be essentially in a single orbit which is that of 
the progenitor, and that orbit happens to be preferentially aligned 
for its transit detection. That is, it is not excluded that we have ob¬ 
served all the clumps present in the system. While a quick look at 
Eigure 11 suggests that the lack of infrared excess might still be 
problematic for the closest pericenters (noting that atot also needs 
to be increased by ^/2 due to the higher transverse velocity at peri¬ 
center in Equation 5), in fact that is not necessarily the case. Rather, 
in that figure we assumed that the clumps were present at the given 
distance at all times, whereas the clumps in the comet-like group 
scenario were at much larger separation from the star at the time 
of the WISE observations. The total mass of the fragmented body 
was considered in Section 4, but since the clumps can be closer to 
the star in this scenario, and are moving faster than circular Keple- 
rian velocity, a better minimum mass estimate for clumps seen at a 


pericenter of 0.1 AU is ~ 3 x 10^® g. Again, the size distribution 
and any material not contributing to the observed dips will increase 
this minimum mass, perhaps by a factor of 100, leading to a more 
realistic parent body mass of 3 x 10^^ g , consistent with a rocky 
body ~ 100km in diameter. 

It remains to be shown that this model can explain the more 
detailed structure of the light-curves. Some potential positives are 
that the clustered nature of the dips could be explained by sub¬ 
sequent fragmentation of a large fragment from an earlier break¬ 
up. The smaller dips could also potentially be explained by smaller 
fragments which may also be expected to receive larger kicks dur¬ 
ing fragmentation. However, the structure of individual clumps may 
be problematic. Eor example, a fairly generic prediction of transits 
of comet-like bodies may be that their light-curves show signs of 
their tails. The light-curve expected for a typical event then has 
a relatively fast ingress as the head of the comet passes in front 
of the star, but a slower egress as the tail passes (e.g. Lecavelier 
Des Etangs et al. 1999; Rappaport et al. 2012). However, the D800 
event shows the opposite (see panel ‘c’ in Figure 1). Possible res¬ 
olutions of this issue are that the D800 comet fragment received 
a large kick with an orientation that sheared it out in such a way 
to form a “forward tail”. Such forward comet tails produced by 
the fragments being kicked toward the star have been studied in 
the literature, but require the grains in the tail to be large enough 
to overcome the effects of radiation pressure (Sanchis-Ojeda et al. 
2015). Alternatively, this event could be comprised of two dips su¬ 
perimposed to have the appearance of a forward tail. While several 
issues remain to be explored, of the scenarios considered we con¬ 
clude that a cometary-like group of bodies seems most consistent 
with the data at hand. 


6 SUMMARY AND CONCLUSIONS 

In this paper, we have shown that KIC 8462852 is a unique source 
in the Kepler field. This otherwise seemingly normal F star un¬ 
dergoes erratic and completely unpredictable dips in flux ranging 
from < 1% to more than 20%. Most of the approximately 7 dips 
observed before D1500 have fairly smooth, but unexplained, dip 
profiles that are each several days long. The D1500 sequence lasts 
continuously for at least 80 days, but the majority of that time is 
spent with the flux depressed by less than ~2%. 

We have conducted numerous follow-up investigations of the 
star and its environment, including spectroscopy, adaptive optics 
imaging, construction of a spectral energy distribution, generation 
of a Fourier transform and a sonogram using the Kepler time series, 
and examination of ground-based photometry. Our analysis charac¬ 
terizes the object as both remarkable (e.g., the “dipping” events in 
the Kepler light curve) and unremarkable (ground-based data reveal 
no deviation from a normal F-type star) at the same time. 

An extensive set of scenarios has been presented to explain 
the occurrence of the dips, most of which are unsuccessful in ex¬ 
plaining the observations in their entirety. Of these, the scenarios 
invoking intrinsic variability, such as the Be star framework, were 
deemed unlikely, but they are not entirely ruled out as a plausible 
option to explain the dips. However, we pointed out that the rel¬ 
atively low Tefj and lack of Ha emission and IR excess in KIC 
8462852 are not suggestive of Be-star activity. 

A broad range of scenarios for the dipping behavior that in¬ 
volve occultation by circumstellar dust clumps was considered. 
Among these, we find that the break-up of one or more massive 
exocomets (or planetesimals on comet-like orbits) provides the 
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most compelling explanation consistent with the data in hand. The 
required mass of the original body may have been in excess of 
3 X 10^^ grams (only ~0.3% the mass of Ceres, and perhaps ~100 
km in diameter). 

We can envision a scenario in which a barrage of bodies, such 
as described above, could be triggered by the passage of a field star 
through the system. And, in fact, as discussed in Sect. 2.3, there is a 
small star nearby (~ 1000 AU; Section 2.3) which, if moving near 
to KIC 8462852, but not bound to it, could trigger such a barrage 
into the vicinity of the host star. On the other hand, if the companion 
star is bound, it could be pumping up comet eccentricities through 
the Kozai mechanism. Measuring the motion/orbit of the compan¬ 
ion star with respect to KIC 8462852 would be telling as to whether 
or not they are physically associated, and we could then be better 
able to make assessments about the timescale and repeatability of 
comet showers based on bound or unbound star-comet perturbing 
models. 

Continuing observations of KIC 8462852 should aid in unrav¬ 
eling the peculiar dips in its light curve. First and foremost, long¬ 
term photometric monitoring is imperative in order to catch future 
dipping events. It will be helpful to know whether such observa¬ 
tions reveal continued, possibly periodic dips, or no further dips. If 
the dips continue, it will be important to search for a clear periodic¬ 
ity, and to look for changes in depth or shape. To completely solid¬ 
ify the hypothesis that the dips are due to dust, observations should 
study the wavelength dependence of the obscuration soon after a 
new dip is discovered. In the case of a family of giant comet-like 
bodies there presumably should be at least a few events similar to 
those seen with Kepler over the next decade. However, if the comet¬ 
like objects actually populate a very long eccentric orbit (i.e., that 
of the original planetesimal), the material may be spread out around 
that orbit, and future dippings events could continue to appear over 
hundreds of years. 

Several of the proposed scenarios are ruled out by the lack of 
observed IR excess (Section 2.4), but the comet/planetesimal frag¬ 
ments scenario has the least stringent IR constraints. In the comet 
scenario, the level of emission could vary quite rapidly in the near- 
IR as clumps pass through pericenter (close to the time they are 
transiting) and are shedding new material. If the system is currently 
in the aftermath of a giant impact, there could be a semi-steady in¬ 
crease in IR flux over years/decades. The WISE observations were 
made in Q5, and assuming that an impact occurred in Q8 (D800, 
Section 5.2), detecting the IR emission from such an impact is still 
a possibility in the future. The only Spitzer IRAC observation of 
KIC 8462852, taken in January 2015, showed a marginal, but below 
3 -(t, excess at 4.5^m, disfavoring the impact scenario (Marengo 
et al. 2015). Continued monitoring in the IR will allow us to firmly 
distinguish between the giant-impact and cometary-group scenar¬ 
ios. 

In summary, it will require some observational skill and pa¬ 
tience to find the next dipping event from this object using ground- 
based observations. As we pointed out, the source spent a rather 
small fraction of its time during the 4-year Kepler mission with dips 
of greater than 2%. Nonetheless, the key to unraveling the mysteri¬ 
ous dips will require such observations. 
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